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Abstract 

Existing theories of migration either focus on micro- or macroscopic behavior of populations; that is, either the 
average behavior of entire population is modeled directly, or decisions of individuals are modeled directly. In this 
work, we seek to bridge these two perspectives by modeling individual agents decisions to migrate while accounting 
for the social network structure that binds individuals into a population. Pecuniary considerations combined with 
the decisions of peers are the primary elements of the model, being the main driving forces of migration. People 
of the home country are modeled as nodes on a small-world network. A dichotomous state is associated with each 
node, indicating whether it emigrates to the destination country or it stays in the home country. We characterize 
the emigration rate in terms of the relative welfare and population of the home and destination countries. The time 
evolution and the steady-state fraction of emigrants are also derived. 

I. Introduction 

There exists no prevalent and coherent theory of international migration yet. Examples of notable observations 
in this area include existence of ethnic clustering, changes in migration patterns (such as Europe, in the second 
half of the 20th century) and occasional overflows of immigration volumes (such as Mexico-US) p4| . This work 
considers the structural processes active in emigrant-emanating and immigrant-absorbing countries, and the micro 
processes motivating individuals and enhancing mobility. 

First efforts to lay a theoretical foundation to study emigration date back to 1889 (Tj. Presently, various approaches 
are segmented by difference in perspective, level of analysis and theoretical tools. Neoclassical economics on the 
macro level reduces migration to wage differentials, whose elimination wipes out the migration flow. It makes labor 
markets the driving mechanisms of international migration, which should be focused on to control the in/out-flow 
of labor ||2]-||5). This approach is accompanied by equivocal empirical investigations |6j, J7J. 

On the micro level, studies are still mainly based on monetary features, where agents optimize productivity given 
their skills and necessary costs. The expected net return at the outset is modeled as follows ||8")-[flT| 

e- rt dt-C(0) (1) 

where P\(t) is the probability of evading deportation from the destination country, Pzit) is the employment 
probability at the destination, Y^t) is the average wage there, P^t) is the employment probability in the home 
country, Yo(t) is the average earning in the home country, r is the rate of inflation, n is the time horizon, and C(0) 
is the value of total migration costs at the outset. Migration goes on as long as the expected net return is positive. 
This model is supported by empirical evidence ]T2)-p4). 

Both of the approaches described above fail to capture the topology of the social network people live in. Another 
potent factor is the stock of immigrants already present at the destination, providing beneficial network externalities 
(the influence of social peers), performing as conduits of information and support, facilitating adjustment to the 
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new place fl4) , fT6)-|[T9). This feature is backed by ample empirical evidence |20|-p5). Network externality and 
cascading (or herding) effect are two inevitable elements of the network perspective. 

In this contribution we approach migration processes from the perspective of network science. We study the 
effect of social influence and wage differentials on the fraction of the population who emigrate. We also introduce a 
bias towards the home land, to model the inherent patriotism people usually have, and find its effect on emigration 
rates. 

II. The Model 

To model social interactions within the population of a country, we employ a simple small-world network 
model p6| . At first, nodes are configured symmetrically on a circle, and each node is connected to every node 
within the proximity of b steps away from it, giving them an initial degree of 2b. Then each pair with distance 
more than b is connected with probability p. The expected total number of links is then Nb + Ar ( jv ~ 1 ~ 2b ) p an( j t ne 
average degree of nodes will be 26 + Q(Np). To ensure that the average degree is finite in the limit as N — > oo it 
is required that p ~ . 

Each node x is endowed with a state s x (t) G {+1,-1}, indicating whether the node stays in its home country 
(+1) or emigrates to the destination country (—1). The states evolve over time according to dynamics which depend 
on factors discussed next. At each time step, a node observes the state of the adjacent nodes and finds the fraction 
of those that are +1 or —1. Denote the fraction of the entire population who, at t, will stay home (with s x (t) = +1) 
by p(t); emigrants (with s x (t) = — 1) are the other 1 — p(t) portion. Each node also observes the expected wage in 
the home country and the foreign country. Here we consider a simplified case where the wage of a country equals 
a characteristic indicator of its national wealth, denoted by G, divided by its population. We denote the population 
of the home and destination country by Nh and N d respectively, and define a = j^. 

A. Migration Without Patriotic Bias 

We first develop a model in which agents make emotionless, rational decisions solely on social interactions and 
pecuniary considerations (sometimes referred to as homo economicus), without any bias towards their homeland. 
The expected wage per capita for the home country is 



and for the destination country it equals 



W - [t) = N d + N G h {l- P (t)) (3) 



If one was to only consider these two criteria (wages in the home and destination countries), the probabilities 
corresponding to leaving and staying would become 



r+(«) 



r_(i) 



G 

w+(t) _ iv^(t) _ G h (a + 1 - p(t)) 



w + (*) + »-(*) W^T) + SK+lfegB G ^ a + 1 - '(*)) + Gdp{t) 

W-{t) _ N d +N h (l-p(t)) _ G d p{t) 



(4) 



We also wish to model social influence. Let M x denote the set of neighbors of node x, and let z x = \Af x \ denote 
the degree of node x. We model social influence in a manner akin to the conventional voter model J28| . Each node 
adjusts its decision (in expectation) according to the fraction of neighbors who are staying (i.e., those neighbors 
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y € M x with s y (t) = +1) and those who are leaving (i.e., s y (t) = — 1). The probabilities it assigns to the actions 
of staying or going are proportional to these fractions. Combining this with the wage considerations, for node x 
we get the following decision probabilities 



P{s x {t + At) = +1} = 



+ (l-0)r+(t) 



(5) 



P{s x {t + At) = -1} = 



1 - 



+ (l-0)r_(t) 



where g (0, 1) is a parameter quantifying the weight given to network effects vs. monetary considerations, and 
5 (a, b) is equal to unity when the arguments match and is zero otherwise. 

We would like to characterize the dynamics of p(t). To do this, we study how the average state evolves with 
time. The expected value of the state for node x in the next timestep is 

E{s x (t + At)} = {+l)P{s x (t + At) = +1} + (-l)P{s x {t + At) = -1} = 2P{s x (t + At) = +1} - 1. (6) 

Using <[3j we get 



E{s x (t + At)} = 26 



2M - 0) G h (a + l-p(t)) _ 
[ >G h {a + l-p{t)) + G d p{t) ■ 



(7) 



Now we need to derive an expression for the first term. Let A denote the adjacency matrix of the network. Then 
we have 

EyeJV, S ( s y: 1) _ A xy S(Sy, 1) 



E 

y=£x 



(8) 



In the so-called annealed approximation ||27j, one studies quantities averaged over the ensemble of all graph real- 
izations (where the average of some quantity -0 denoted by (tp)) , and approximates j4 - x "^ s "- 1 ) by (^r L )(5{s yi 1)). 
All AijS are random variables with binary distribution, and the average is over all configurations. We need to 
consider two different cases, depending on whether or not x and y are initially within b hops of each other in the 
small-world model, before long-range connections are added. If y is one of the 2b initial neighbors of x, then 



*xy 



N 

E 



P{z x = k + 2b} 
k + 2b 



E 



ft ) p Hi-pY 

k + 2b 



def 



(9) 

.. , _ /.' -+ n> 

k=0 k=0 

where N' =N — 2b—\ is the maximum number of links which may be added, per node, in the small-world model. 
If y is not initially linked to x then we have 



.4 



■nj 



_ P{z x = k + 2b + l} ^ (V)P»(1-/- 
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= 0. 



(10) 
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As noted in \ 21\, one can express uj as the following integral 



and similarly 



Integrating by part gives 



from which it is clear that 
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Combining this with (j7]i we get 

2p(t + At) - 1 = 20p(i) -1 + 2(1 
which leads to the following differential equation for p 

p=(0-l)p+(l- 



2bcu + n(N -2b- 1) = 1. 

G h (a + l-p(t)) 



'G h (a + 1- p(t)) + G d p{t) 
G h (a + l-p) 



G h (a + l- p) + G dP 
Note that this differential equation is separable, and so its solution is expressed as 

G h (a + l-p) 



rp(t) 

(l-9)t= / dp 

Jp{0) 



-P 



G h (a + 1 - p) + G dP 



(11) 
(12) 

(13) 
(14) 

(15) 
(16) 

(17) 

(18) 
(19) 
(20) 
(21) 

(22) 
(23) 
(24) 
(25) 
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The closed form expression for this integral is given by 



1 



In 



(1 + a)G h (l - pjt)) - p(t)(l - p{t))G h - G d p{tf 
—Gd 

G h (a+2-2p(t))+2G d p(t) 

/G^ 1 /G h (l+a) 2 +4(l+a)G d -G fe (l+2a) 



a\/Gh tanh 1 



aJGh tanh 



G h a+2G d 



/G^^G h (l+a) 2 +4(l+a)G d -G h (l+2a) 



y/G h {l + a) 2 + 4(1 + a)G d - G h {\ + 2a) 



(26) 



which is not very informative. Thus we seek the steady-state solution which can be found by setting the time 
derivative equal to zero, which gives 

def .. ... V(a + 2)2 + 4G h (G d ~G h )(l + aj~(a + 2) 

Poo = lim pit) = — — . (27) 

t^oo 2(Grf — &h) 

We validate the expressions derived above via simulation. Figure [T] depicts the equilibrium fraction of emigrants 
as a function of a for different values of G d when Gh = 1- Figure [2] shows the expected probability of emigration 
as a function of the fraction of those who already migrated. 




4 5 6 7 

Population ratio a 



Fig. 1: The steady-state fraction of the population who decide to leave (1 — Poo), as a function of a, for different 
values of G d , when Gh is normalized to unity, and 8 = 0.7. There are N = 1000 nodes, with 6 = 8 and p = 4/N, the 
steady-state is acquired after 2000 timesteps. Solid lines and bullets represent theoretical prediction and simulation 
results respectively. 



B. Migration With Patriotic Bias 

Note that so far there is no intrinsic tendency towards "home", homo-economicus-type agents are in action and 
high probabilities are not surprising. For instance, this helps explain why the fraction of nodes emigrating can reach 
70% in Figure [T] 

To be more realistic, one must account for the disposition against emigrating due to patriotism. We account 
for this in a Bayesian manner. Let 1 — 7 denote the prior probability of staying at home regardless of social and 



July 26, 2012 



DRAFT 



6 




Fig. 2: Emigration probability P{s x (t + At) = — 1} as a function of the fraction of those already emigrated, 
1 — p(t), for different values of Gd, when Gh = 1. There are N = 1000 nodes, with b = 8 and p = A/N . 



economic conditions. Incorporating this factor into the model described in the previous section leads to the expected 
state dynamics 



E{s x (t + At)} = 7 



{+l)P{s x (t + At) = +1} + {-l)P{s x {t + At) = -1} 



+ (l-7)(+l) 



2P{s x {t + At) = +!}-! 



+ 1- 7- 



The equation of motion becomes 



p = (70-l)p + 7(l-. 



G h {a + \-p) 



G h (a + 1- p) + G dP 
The solution to this differential equation can be found by taking the integral 

dp(t) 



I-7. 



t 



(28) 



(29) 



(30) 



(31) 



(i - 7) + prow - 1) + 7(1 - o) ^m^mm 

which is straightforward but voluminous. One can simplify the above expression when the fraction of emigrants is 
very smalQ by setting p(t) = 1 — e, and taking a Taylor series expansion about e = up to second order terms. 
Defining 

'« s 7(l-«)( 5n ?fei) 

^l-70 + t^(l + ^g^) (32) 

= 7(l-g)Gh(Gh-Gd)(a-l) 



'According to the CIA World Fact Book, the vast majority of countries have migration rates below 1%. 
https://www.cia.gov/library/publications/the-world-factbook/rankorder/21 12rank.html 
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we have 



t 



1 



C-7^ 2 +(^-f )< :i Ol,') 



Alternatively, the equilibrium emigration fraction is found by solving 

G h (a + 1- P (t)) 



( 7 fl-l)p(t) +7 (l- 



whose solution is 



with 



Poo = lim p(t) = 

t— >oo 



G h (a+l-p{t)) + G d p(t) 
-B - V-B 2 - 4AC* 



1-7 = 



2/1 



A=(l- 7 0)(G, i -G d ) 

S = (a + l)G h ( 7 6 - 1) + (1 - 7)(G d - - 7 G h (l - 0) 
C=G h (l + a)(l-70). 



(33) 



(34) 



(35) 



(36) 



Note that the sign of poo depends on the sign of A; if Gd > Gh, then p^ > and there is emigration to the 
destination country, as one would expect. 

The equilibrium fraction of emigrants is depicted in Figures [3] and [4] and the emigration probability with respect 
to the fraction of those already done so is given in figure |5]), which replicates the one due to network externalities 
and herding effect in fT7) . 




Fig. 3: The fraction of the population who decide to leave, as a function of a, for different values of Gd, when Gh 
is normalized to unity, with 9 = 0.7, 7 = 0.2, N = 1000, 6 = 8 and p = 4/JV. Solid lines and bullets represent 
theoretical prediction and simulation results respectively. 



Note that in the calculations, the role of p (the connection probability the small-world model) cancelled out. In the 
simulations no effect on the emigration fractions was observed by changing p. This peculiarity of the small-world 
network coincides with the findings in [27 1, where the average magnetization is preserved and the correlation length 
of the finite graph reaches unity in the steady state, irrespective of the value of p. When there is perfect symmetry 
between all nodes, local densities vanish in the long time limit, and the equilibrium fraction of emigrants (or in the 
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Fig. 4: The fraction of the population who decide to leave, as a function of a, for different values of 9, when G/j 
is normalized to unity, with Gd — 2, 7 = 0.2, N = 1000, 6 = 8 and p = 4/N. Solid lines and bullets represent 
theoretical prediction and simulation results respectively. 



voter model as in (27], average magnetization) only depends on the total fraction of the emigrants. This suggests 
that, in order to further incorporate the structural properties of the society, one must introduce asymmetry between 
nodes by construction, that is, enforce hierarchy of influence or stratify the underlying network. 

Also note that, in the case of 7 = 9 = 1 we have the conventional voter model, which only has two absorbing 
states; either all nodes end up with s = +1 or with s = — 1. Figure [6] represents the probability of al the population 
staying in the home country, as a function of those who decided to do so at the outset. Note that the results are 
the same for different values of p. 



Figures 7a and 7b show the effect of 9 on the fraction of emigrants. The network has N — 80 nodes with 
connection probability p = 2/N, and all other parameters rather than 9 are the same for the two graphs. 

III. Summary 

Previous models proposed for the emigration process do not take into account the social network structure of 
the population. In this paper we have modelled emigration decisions taken by agents interacting over a small-world 
network. The first model was merely utility based. The fraction of emigrants were found in terms of the difference in 
average wage in the home and destination countries, the population ratio of the population of the two countries, and 
the relative importance that nodes give to their personal utility calculations and to imitating the behavior of others. 
The second model embodied a default disposition against leaving the homeland, modelling patriotism. Theoretical 
predictions were accompanied by simulations. We also found that as long as the nodes are indistinguishable, i.e. 
their position and characteristics in the network are identical, then the average emigration rates are independent of 
the connection probability in the small-world network. One must introduce hierarchy or other structural asymmetries 
in order to study the dependence of the collective behavior on the structural properties of the underlying network. 
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Fig. 5: Emigration probability , as a function of 1 — pit), the fraction of those already emigrated, for different 
values of Gd, when Gh is normalized to unity, with 7 = 0.2. 
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is set to 0.2, and for the one on the right, it is set to 0.8. 
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